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Abstract 

We review self-consistent spectral methods for nuclear matter calculations. The 
in-medium T— matrix approach is conserving and thermodynamically consistent. It 
gives both the global and the single-particle properties the system. The T— matrix 
approximation allows to address the pairing phenomenon in cold nuclear matter. A 
generalization of nuclear matter calculations to the superfiuid phase is discussed and 
numerical results are presented for this case. The linear response of a correlated sys- 
tem going beyond the Hartree-Fock-I- Random-Phase- Approximation (RPA) scheme 
is studied. The polarization is obtained by solving a consistent Bethe-Salpeter (BS) 
equation for the coupling of dressed nucleons to an external field. We find that 
multipair contributions are important for the spin(isospin) response when the inter- 
action is spin(isospin) dependent. 



1 Introduction 



Nuclear matter is an infinite system of strongly interacting fermions. Low energy nucleon- 
nucleon interactions in vacuum are well known from scattering experiments. The char- 
acteristic features of the nuclear force are the presence of a strongly repulsive core at 
small distances, attraction at moderate distances, and the appearance of a tensor com- 
ponent. The properties of nucleons in a strongly interacting medium are modified. This 
is taken into account by the dressing of nucleon propagators by a self-consistently cal- 
culated self-energy. Nonzero imaginary part of the self-energy implies a finite lifetime of 
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quasiparticles and gives them nontrivial spectral properties; nucleons propagate off the 
energy shell. The description of nuclear matter in the language of such dressed nucleons 
with a broad spectral function has been the subject of intensive studies in the last years. 

2 In medium T matrix 

The description of properties and excitations of a many-body system can be performed 
in the language of finite-temperature Green's functions. The Green's functions' approach 
is most easily formulated in the imaginary time formalism P, suitable for formal presen- 
tations and perturbative calculations, but the real-time formalism [2] is better adapted 
for numerical calculations. In medium propagators (Green's functions) are dressed by a 
self-energy term. The physical approximation enters in the choice of a suitable expression 
for the self-energy used in the calculation. Analytical properties of the Green's functions 
are guaranteed by the dispersion relation between the imaginary and real parts of the 
self-energy and by the use of the Dyson equation 

G- 1 = Go 1 - S , (1) 

where G and G represent the dressed and vacuum propagators respectively and £ is the 
self-energy. 

The choice of a specific form of the self-energy should be adapted to the physical system 
under consideration, but a general procedure to derive nonperturbative approximations 
in many-body systems has been proposed The self-energy is written as a functional 
derivative of a two-particle irreducible generating functional S = 

Due to the presence of a short range repulsive core, ladder diagrams in the free nucleon- 
nucleon interaction V must be resummed. The resulting correlated part of the two-particle 
Green's function is called the in-medium T— matrix jl] (Fig. ^) 

T = V + VGGT . (2) 

Although the ladder diagrams for the T matrix have the same form as for the Brueckner 
G— matrix the resulting expressions are different. The T— matrix equation is defined for 
in-medium Green's functions. At zero temperature the two-nucleon propagator in Fig. ^ 
represents the propagation a particle-particle pair (excitations above the Fermi energy) or 
a hole-hole pair (excitations below the Fermi energy), whereas in the Brueckner scheme 
the blocking operator forces the two nucleons in the ladder to be of the particle-particle 
type. The second, even more important, difference comes at the level of the self-energy, 
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Figure 1: Definition of the T— matrix as a sum ladder diagrams in the interaction. 




Figure 2: The self-energy in the T— matrix approximation 

which in the $— derivable T— matrix approximation must be taken self-consistently in the 
form (Fig 

S = TrTG . (3) 

The self-consistency means that the nucleon propagators at all stages of the approxima- 
tion are dressed by this self-energy. It requires an iterative procedure for the solution of 
the coupled equations (J2J), ©, and (JTJ). The T matrix obtained in a given iteration is 
used to generated the self-energy and the dressed propagators. The dressed propagators 
are then used to calculate the in-medium T matrix in the next iteration. 

The numerical calculations are very demanding due to the fact that the nucleons are 
dressed by nontrivial spectral functions coming from the imaginary part of the self-energy 
(FigEl); additional integrations over the energy of off-shell nucleons appear. First results 
on the self-consistent in-medium T matrix appeared only in the recent years El • in 
El the spectral function of the dressed nucleon was parameterized by three Gaussians 
and in [HI E] using a three poles approximation for the in-medium propagator. The first 
solution using full spectral functions for the dressed propagators in nuclear matter was 
obtained numerically in ref. [0] at finite temperature, for a simple nuclear interaction. 
This approach, using numerical algorithms significantly speeding the calculations, was 
then generalized to realistic interactions and zero temperature [TU1 ITT] . The solution of 
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Figure 3: The nucleoli spectral function obtained in the self-consistent T— matrix approx- 
imation for nuclear matter. 

the T— matrix scheme at finite temperature was also obtained in |12j . 

At low temperatures two difficulties show up. The first one is related to the transition 
to superfluidity in cold nuclear matter. This effect can be taken into account quite nat- 
urally within suitably generalized T— matrix approaches. This is the subject of the next 
section. We note however that most of the existing T— matrix calculations are done for 
the normal phase of the nuclear matter also at zero temperature. The second difficulty 
with low temperature nuclear matter is technical and is related to the appearance of a 
well defined quasi-particle peak in the spectral function of the in-medium nucleon. Simple 
discretization algorithms for the energy integration break down in that case. The spectral 
function must be separated into background and quasi-particle contributions for momenta 
close to Pf- At the time of this writing, there is only one numerical implementation of 
this procedure for the nuclear matter T matrix (see JU] for details). 

First calculations Ej have shown that the self-consistency for the self-energies 
of nucleons is very important. The effective scattering is reduced in the self-consistent 
calculation when compared to an approximation neglecting the imaginary part of the self- 
energy Also the value of the critical temperature for the superfluid phase transition 
goes down when using fully dressed propagators in the Thouless criterion for superfluidity. 

Further developments have addressed the binding energy in the T— matrix approxima- 
tion flj EE El OS] . The role of the correlations, high momentum states, and low energy 
tails in the spectral function on binding has been discussed [T3]- Formally, one should 
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Figure 4: The binding energy and the Fermi energy (upper panel) and two (formally equiv- 
alent) expressions for the pressure (lower panel) as function of density in the T— matrix 
and the Brueckner-Hartree-Fock calculations. 



not expect better results for the binding energy in the T— matrix approximation than in 
the standard G— matrix approach with higher order terms in the hole line-expansion [T7j . 
However, the authors of ref. ^H] claim that the T— matrix calculation, which does not 
take into account ring diagrams contribution, is better adapted to describe nuclear mat- 
ter in finite nuclei. An important observation is made in ref. ^3], where the role of the 
thermodynamical consistency of the T— matrix approximation is discussed. As expected 
from the $— derivability of the scheme, the T— matrix calculation gives consistent results 
for the global quantities in the system and its single particle properties (Fig. |3J). In par- 
ticular the T— matrix results fulfill the Hugenholz-Van Hove theorem and the Luttinger 
relation for the value of the Fermi momentum. In practice, such a consistency of the 
single particle energy can be obtained within the Brueckner theory only approximately 
by taking rearrangement terms to a given order. 

The single-particle properties obtained in the self-consistent iteration of the T— matrix 
scheme represent a reliable estimate of the optical potential and of the single-particle width 
in medium [IjQj. Nuclear matter in the non-superfluid phase behaves as a Fermi liquid, 
with the standard phase-space scaling of the scattering width with the distance to the 
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Fermi energy and with the temperature • The spectral function in finite nuclei has been 
estimated taking into account the local density and the neutron-proton asymmetry ^H] 
and compared to experimental data. Thermodynamic properties of the asymmetric 
and pure neutron nuclear matter [20J have been discussed. We note that the T— matrix 
calculation can be performed very easily at finite temperature [31101112] (it is even simpler 
than at zero temperature). 

3 Pairing 
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Figure 5: The spectral function in the superfluid phase, including the the diagonal self- 
energy (solid line) and both the diagonal and anomalous self-energies (dashed line). 

A general property of fermions interacting through an attractive potential is the tran- 
sition to a superfluid state at low temperatures pT| . This phenomenon is usually not 
taken into account in nuclear matter studies, because the value of the superfluid gap, as 
extracted from the properties of finite nuclei, is expected to be small. On the other hand 
naive estimates of the neutron-proton pairing gap yield a value of several MeVs j22] • This 
value, obtained from a mean-field gap equation, is modified by many-body effects. 

A fundamental approach in the study of superfluid nuclear matter should aim at the 
description of the cold nuclear matter system dealing with a strong repulsive core in the 
interaction and at the same time with the formation of the superfluid state. A generaliza- 
tion of nuclear matter calculations must include the usual ladder diagram resummation 
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in the normal phase. At the same time it should describe the formation of the superfTuid 
long-range order in the two-particle correlations. Such correlations appear as a singularity 
of the T matrix for T = T c at the energy of twice the fermionic chemical potential. This 
corresponds to the formation of a two-fermion bound state. At temperatures below T c 
such fermionic (Cooper) pairs have a binding energy allowing for the condensation and 
for the formation of a superfluid order parameter. 
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Figure 6: Inverse of the T— matrix in the pairing (deuteron) channel in the superfluid 
phase. The usual T matrix (dashed line) and the generalized T matrix with off-diagonal 
components (solid line) are shown |2Sj- The generalized T matrix has a singularity at the 
Fermi energy for zero total momentum of the pair. 

In ref. a first attempt to generalize nuclear matter calculation to the superfluid 
phase has been discussed. At small temperatures a superfluid state forms leading to 
nonzero expectations of the off-diagonal (anomalous) Green's functions and self-energies 
[25] . The off-diagonal self-energy is obtained from a gap equation but with dressed prop- 
agators. The T— matrix in the superfluid must be constructed in a way to conserve the 
singularity at twice the chemical potential also for T < T c . Such a scheme was considered 
in [21], giving the first calculation of the superfluid nuclear matter problem, including 
both the ladder resummation of the nuclear interaction and the superfluid properties ob- 
tained with a gap equation modified by many-body effects. The nucleon propagator in 
the superfluid is dressed by a diagonal self-energy S coming from a modified T— matrix 
approximation and an off-diagonal anomalous self-energy A; it has two poles on both 
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sides the Fermi energy (Fig. |5J). The most interesting result of this first, and up to now 
unique, study is the observation of a strong reduction of the superfluid gap in the cor- 
related nuclear matter as compared to results from a mean-field gap equation [21]. This 
effect was analyzed in details [26J and can be explained by an effective modification of 
the density of states at the Fermi energy. The reduction of the superfluid energy gap 
is especially important for neutron-proton pairing close to the gap closure in symmetric 
nuclear matter [77] . 

Possible generalizations of the T— matrix approximation to the superfluid phase have 
been discussed by Haussmann [2E]- The simplest scheme, which is derivable and 
includes the required ingredients, i.e. the ladder diagrams resummation in the diagonal 
self-energy and the gap equation for the off-diagonal self-energy, requires the introduction 
of an additional T— matrix describing two-particle correlations for anomalous propagators. 
The generalized T— matrix has a singularity at twice the Fermi momentum for T < T c 
(Fig. IHl) and the imaginary part of the self-energy has a gap around the Fermi energy, 
corresponding to an energy gap for possible excitations in a scattering [231- This last, 
important property of the approximation does not hold for the simple scheme discussed 
in |23|. Explicit calculations show that corrections to the gap equation coming from ladder 
diagrams in the off-diagonal self-energy are small [23] . These small corrections make the 
superfluid gap energy dependent. 

4 Linear response with dressed vertices 

Processes occurring in the dense nuclear matter are modified by medium effects 
This happens for neutrino rates in neutron stars and for particle or photon emission in 
hot nuclear matter. In particular, the effect of the off-shell propagation of nucleons in the 
medium is important for the subthreshold particle production in heavy ion collisions. The 
role of correlations for the neutrino emission could be especially important for processes 
in hot stars. 

If the interaction with an external perturbation is small the dynamics of the system 
in the external field can be described in terms of linear response functions. The response 
functions incorporate all the correlations due to the self-interaction of particles in the 
system. For normal Fermi liquids the response function to long-wavelength perturba- 
tions can be calculated within the Fermi liquid theory j3Uj. However, if one is interested 
in nuclear systems at higher temperatures or for perturbations with large momentum a 
different approach is required. The Fermi liquid theory does not allow for multipair exci- 
tations of the system in the external potential. Such multipair excitations are important 
for higher energies and momenta of the external perturbation or in the presence of tensor 
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Figure 7: Diagrams for the self-energy. The first two diagrams are the Hartree-Fock 
contribution for the Gogny interaction. The last diagram is the contribution of the residual 
interaction in the second order. 

interactions |31j . 

We take the interactions between the nucleons as a sum of a mean-field interaction that 
is based on the Gogny parameterization (221 an d a residual interaction (scalar or isospin 
dependent). The isospin dependent residual interaction is obtained from the scalar one 
multiplying by the factor |(1 + tit 2 ). The self-energy is taken in the second direct Born 
approximation for the residual interaction (Fig. |7J). The residual interaction induces a 
finite width to nucleon excitations in the medium. Such a dressing of nucleons is expected 
in any approach going beyond the simple mean-field, e.g. the T— matrix approximation 
discussed previously. 

When the description of the correlated system goes beyond the mean-field approxima- 
tion the difficulty involved in a consistent calculation of the response function is severely 
increased [331 BH • A naive calculation of the polarization bubble using dressed propaga- 
tors 

IT = TrT G ph T , (4) 

where G p h is the particle-hole propagator with dressed nucleons and r is the free vertex 
for the coupling of the nucleon to an external field, is a very bad estimate for the response 
function [33]. In particular, it severely violates the tu-sum rule [3Sj . 

A general recipe for calculating the in-medium coupling of the external potential to 
dressed nucleons is known [37j. The in- medium vertex describing the coupling of an 
external perturbation to nucleons is given by the solution of the BS equation (Fig. [Sj) 

T = T + TrKG ph T , (5) 

where K denotes the particle-hole irreducible kernel. The kernel K of the BS equation 
should be taken consistently with the chosen expression for the self-energy, it is given by 
the functional derivative of the self-energy with respect to the dressed Green's function 
K = 5T>/5G [37|. Using the dressed vertex obtained as a solution of the BS equation the 
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Figure 8: The Bethe-Salpeter equation for the dressed vertex. The particle-hole irre- 
ducible kernel K is denoted by the box and the fat and the small dots denote respectively 
the dressed and the bare vertices for the coupling of the external field to the nucleon. 
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Figure 9: The polarization function expressed using the dressed vertex for the coupling 
of the external field to the dressed nucleon. 

response function in the correlated medium can be obtained from the diagram in Fig. |§] 

ri(ST) = TrT° ST ^G ph T^ S T) , (6) 

T(sT) is t ne in- medium vertex for the coupling in a given spin-isospin (ST) channel. The 
exact form of the BS equation in the real-time formalism (for different types of vertices 
T(5T)) can be found in [313 EB]- The solution of the BS integral equation is obtained by 
iteration. It is a serious numerical task involving the calculation of two-loop diagrams with 
dressed propagators and broken rotational symmetry due to the presence of an external 
field. 

The results for the response functions in different spin-isospin channels jHH] show that 
the calculation using dressed propagators and dressed vertices obtained as solutions of 
the BS equation is close to the RPA approximation. It is expected due to cancellations of 
self-energy and vertex corrections. Such cancellations occur for the exact solution of the 
system as well as within consistent approximation derived from a generating functional 
$. For the scalar residual interaction the w-sum rule takes the simple form 

-/^ Imn (5T)(q^) = P^, (7) 

in all the spin isospin channels (ST). The above sum-rule severely constraints acceptable 
forms of the response functions in the case of a scalar residual interaction and leads to 
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very small multipair contributions in most of the kinematical regions. Some differences 
can occur only when collective modes are present, as discussed latter. 

The situation is very different for a more general form of the residual interaction. The 
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Figure 10: The imaginary part of the polarization in the RPA approximation (dashed- 
dotted line), from the self-consistent calculation with dressed nucleons and vertices (solid 
line) and the naive one-loop polarization with dressed nucleons (dashed line). All results 
are for an isospin dependent interaction |(1 + t 1 t 2 )V, q = 210MeV, and T = 15MeV. 

kernel of the BS equation depends very much on the isospin channel, e.g. in the channel 
ST = 11 there are no vertex corrections from the residual interaction at all. The propaga- 
tors are dressed by a nontrivial self-energy due to the residual interaction but the vertex 
corrections are small or even absent. Therefore, the cancellation between self-energy and 
vertex corrections, observed for a scalar interaction, can no longer be maintained. This 
has its implications for the uj— sum rule which gets modified for responses with T = 1 
[35] . In the case of the isospin dependent residual interaction, in-medium dressed nucleons 
couple in the same way as free nucleons to isovector potentials. In Fig. El the response 
functions obtained with the isospin dependent interaction are compared to the response 
functions from the Fermi liquid theory. For T = channels the correlated response func- 
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tion is similar to the one particle-one hole response function. On the other hand, the 
isovector response is closer to the naive one-loop result, without vertex corrections. 
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Figure 11: The imaginary part of the polarization when a collective mode is present; 
density response (upper panel) and isospin response (lower panel). The results are ob- 
tained in the RPA approximation (dashed-dotted line), and in the full calculation with 
dressed vertices for an isospin dependent residual interaction (solid line) and a scalar one 
(dashed line). For the density response the result is independent of the type of the residual 
interaction. 

To study the role of multipair configurations on the collective modes we increase the 
value of the Landau parameters. Within the Fermi liquid theory a collective excitation at 
zero temperature is a discrete peak in the imaginary part of the response function. The 
state corresponding to a collective excitation cannot couple to incoherent one particle-one 
hole excitations. At finite temperature such a coupling is possible, it can be calculated 
and the finite temperature width of the collective state is usually small. The collective 
state can acquire a finite width (also at zero temperature) due to the coupling to multipair 
configurations [30]. The description of the damping of collective states from such processes 
goes beyond the Fermi liquid theory. 
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The result of the calculation with self-energy and vertex corrections is compared to 
the RPA result in Fig. ^2 The collective mode in the density response has a large 
width in the full calculation. This corresponds to the coupling of the sound mode to 
multipair excitations giving rise to a finite decay width, even at zero temperature (the 
small width of the RPA collective mode is due to the finite temperature T = 15MeV). 
For the isovector response shown in the lower panel of Fig. ^2 the Fermi liquid theory 
predicts the presence of a well defined collective state (dashed-dotted line). For a scalar 
residual interaction the response function in the correlated system shows a collective state 
at similar energy (dashed line). It has a larger width due to the contribution of multipair 
configurations, analogously as in the density response. On the other hand, an isospin 
dependent residual interaction leads to an isovector response with a large imaginary part 
at high energies (Fig. Illljl . For the chosen interaction we observe the extreme scenario of 
the disappearance of the collective mode when the effects of the residual interaction are 
taken into account (solid line in the lower panel of Fig. [TTj) . Generally we expect a whole 
range of behavior depending on the energy of the collective state and on the strength 
of isospin dependent terms in the residual interaction but always the collective state is 
broader than in the Fermi liquid theory, due to the coupling to multipair configurations. 
The same phenomena are expected also for the spin wave collective state in the presence 
of spin dependent residual interactions. 

5 Summary 

We present a new method of calculations for nuclear matter. The method discussed in 
this work is using in-medium Green's functions' formalism and is based on the summa- 
tion of diagrams with a retarded propagation of a pair of fermions. The sum of such 
ladder diagrams is the in- medium scattering matrix (the T-matrix). Self-consistency in 
this approach means that the fermion propagators in the T-matrix diagrams are dressed 
in a nontrivial self-energy. The self-energy itself is obtained from the sum of ladder dia- 
grams. This self-consistency requires the use of full spectral functions depending on the 
momentum and energy of the particle. The propagation off the mass shell, i.e. the full 
dependence of the nucleon propagator on the energy, is a serious difficulty in numerical 
applications. The progress achieved in the last years allows to perform extensive cal- 
culation of the nuclear matter properties in the T— matrix approach at zero and finite 
temperatures. 

Attractive nuclear interactions lead to the formation of a superfluid phase at low tem- 
peratures. Cooper pairs of fermions are formed, analogously as for the electron supercon- 
ductivity in metals. The second important achievement here reported is the generalization 
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of the ladder diagram summation to the superfluid phase. We discuss equations which 
allow for a simultaneous and consistent treatment of the short range nuclear interactions, 
the bound state formation (Cooper pairs) in the T-matrix and the superfluidity in nuclear 
matter. Model calculations have allowed us to estimate the influence of the superfluid- 
ity on standard many-body effects in nuclear matter. At the same time, the influence of 
many-body corrections on the superfluid gap can be assessed. These corrections, reducing 
the value of the gap, can be described by an effective superfluid gap equation, similar to 
the mean-field Bardeen-Cooper-Schrieffer equation with an effective interaction and with 
a renormalized value for the order parameter. 

The self-consistent T-matrix approximation is a thermodynamically consistent ap- 
proximation, a so called conserving approximation. Single-particle properties obtained 
within this approximation are consistent with global quantities describing the system, 
such as the binding energy or the pressure. Unlike other methods, it allows to obtain 
directly single-particle properties : optical potentials, single-particle widths, and spectral 
functions. These observables can be compared to experimental values. 

The calculation of processes occurring in the dense medium, e.g. particle absorption, 
emission, neutrino cross-sections, requires the knowledge of different kinds of linear re- 
sponses of the system to external probes. In a correlated system, described using dressed 
propagators, the calculation of the linear response is a very difficult task. The problem 
lies in the need for a consistent dressing of the vertices corresponding to the coupling 
of dressed nucleons to an external field. The in-medium vertices are obtained from a 
solution of the BS equation, where the kernel of the equation is derived from the same 
generating functional $ as the self-energy. This procedure guarantees the fulfillment of 
conservation laws and sum-rules. It must be contrasted with incomplete procedures with 
only self-energy effects without vertex corrections. The linear response with self-energy 
and vertex corrections beyond the Hartree-Fock+RPA approximation includes the effects 
of multipair excitations. For fermions interacting with a scalar potential such effects are 
small, except for the damping of collective modes, where the multipair configurations are 
important. The situation is very different for spin or isospin dependent interactions. In 
that case the spin or isospin response is very different form the RPA result; multipair 
effects are essential. 

Applications discussed in this talk are restricted to nuclear systems. We note that the 
same methods can be applied to other many-body fermionic systems, high T c supercon- 
ductors, fermions near the Feshbach resonance, the electron gas. 

This work was partly supported by the Polish State Committee for Scientific Research 
Grant No. 2P03B05925. 



14 



References 

[1] T. Matsubara, Prog. Theor. Phys., 14 (1955) 351. 
[2] L. V. Keldysh, Zh. Eksp. Teor. Fiz., 47 (1964) 1515. 
[3] G. Baym, Phys. Rev., 127 (1962) 1392. 

[4] L. Kadanoff, G. Baym, Quantum Statistical Mechanics, Bejamin, New York, 1962. 
[5] W. H. Dickhoff, Phys. Rev., C58 (1998) 2807. 
[6] P. Bozek, Phys. Rev., C59 (1999) 2619. 

[7] W. H. Dickhoff, C. C. Gearhart, E. P. Roth, A. Polls, A. Ramos, Phys. Rev., C60 
(1999) 064319. 

[8] Y. Dewulf, D. Van Neck, M. Waroquier, Phys. Lett., B510 (2001) 89. 
[9] Y. Dewulf, D. Van Neck, M. Waroquier, Phys. Rev., C65 (2002) 054316. 
[10] P. Bozek, Phys. Rev., C65 (2002) 054306. 

[11] P. Bozek, P. Czerski, Acta Phys. Polon., B34 (2003) 2759-2768. 
[12] T. Frick, H. Miither, Phys. Rev., C68 (2003) 034310. 

[13] T. Aim, G. R6pke, A. Schnell, N. H. Kwong, H. S. Kohler, Phys. Rev., C53 (1996) 
2181. 

[14] P. Bozek, P. Czerski, Eur. Phys. J., All (2001) 271. 
[15] P. Bozek, Eur. Phys. J., A15 (2002) 325. 

[16] Y. Dewulf, W. H. Dickhoff, D. Van Neck, E. R. Stoddard, M. Waroquier, Phys. Rev. 
Lett., 90 (2003) 152501. 

[17] H. Q. Song, M. Baldo, G. Giansiracusa, U. Lombardo, Phys. Rev. Lett., 81 (1998) 
1584. 

[18] P. Bozek, Phys. Lett., B586 (2004) 239. 

[19] T. Frick, H. Miither, A. Rios, A. Polls, A. Ramos, Phys. Rev., C71 (2005) 014313. 



15 



[20] P. Bozek, P. Czerski, Phys. Rev., C66 (2002) 027301. 

[21] J. Bardeen, L. N. Cooper, J. R. Schrieffer, Phys. Rev., 108 (1957) 1175. 

[22] B. Vonderfecht, C. Gerhart, W. Dickhoff, A. Polls, A. Ramos, Phys. Lett., B253 
(1991) 1. 

[23] P. Bozek, Phys. Rev., C65 (2002) 034327. 
[24] P. Bozek, Nucl. Phys., A657 (1999) 187. 

[25] J. R. Schrieffer, Theory of superconductivity, W. A. Benjamin, New York, 1964. 

[26] P. Bozek, Phys. Rev., C62 (2000) 054316. 

[27] P. Bozek, Phys. Lett., B551 (2002) 93. 

[28] R. Haussmann, Z. Phys., B91 (1993) 291. 

[29] J. Knoll, D. N. Voskresensky, Annals Phys., 249 (1996) 532. 

[30] D. Pines, P. Nozieres, The Theory of Quantum Liquids Vol. /, Benjamin, New York, 
1966. 

[31] E. Olsson, C. J. Pethick, Phys. Rev., C66 (2002) 065803. 
[32] D. Gogny, R. Padjen, Nucl. Phys., A293 (1977) 365. 
[33] N. Kwong, M. Bonitz, Phys. Rev. Lett., 84 (2000) 1768. 
[34] S. Faleev, M. Stockman, Phys. Rev., B66 (2002) 085318. 
[35] P. Bozek, Phys. Lett., B579 (2004) 309. 

[36] D. Tamme, R. Schepe, K. Henneberger, Phys. Rev. Lett., 83 (1999) 241. 

[37] G. Baym, L. KadanofT, Phys. Rev., 124 (1961) 287. 

[38] P. Bozek, J. Margueron, H. Miither, Ann. Phys., 318 (2005) 245. 



16 



